Spatio-Temporal Change of Land Use/Land Cover and Vegetation Using Multi-MODIS Satellite Data, Western Ethiopia

Land use and land cover (LULC) change and variability are some of the challenges to present-day water resource management. The purpose of this study was to determine LULC and Normalized Difference Vegetation Index (NDVI) fluctuations in western Ethiopia during the last 20 years. The first part of the study used MODIS LULC data for the change analysis, change detection, and spatial and temporal coverage in the study region. In the second part, the study analyzes the NDVI change and its spatial and temporal coverage. In this study, The Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data were applied to determine LULC and NDVI changes over four different periods. Evergreen broadleaf forests, deciduous broadleaf forests, mixed forests, woody savannas, savannas, grasslands, permanent wetlands, croplands, urban and built-up lands, and water bodies are the LULC in the period of analysis. The overall classification accuracy for the classified image from 2001 to 2020 was 85.4% and the overall kappa statistic was 81.2%. The results indicate a substantial increase in woody savannas, deciduous broadleaf, grasslands, permanent wetlands, and mixed forest areas by 119.6%, 57.7% 45.2%, 37%, and 21.3%, respectively, followed by reductions in croplands, water bodies, savannas, and evergreen broadleaf forest by 90.1%, 19.8%, 13.2%, and 4.8%, respectively, for the catchment between 2001 and 2020. The result also showed that the area's vegetation cover increased by 64% from 2001 to 2022. This study could provide valuable information for water resource and environmental management as well as policy and decision-making.


Introduction
Land use and land cover (LULC) change is a key driver of environmental change and is becoming a global concern due to its impact on the local, regional, and global environment [1][2][3][4][5].Understanding how LULC classes afect upcoming precipitation as well as surface dynamics which are ecologically relevant is crucial for sustainable water resource management and ecosystem [6,7].LULC change is exacerbated by human activities such as population growth, urbanization, and deforestation [8].Several studies have identifed these changes as a current challenge to the ecohydrological system of the environment [8][9][10][11][12].
Tere is presently a wide spectrum of studies on LULC, on the database, taking place throughout the world utilizing various technological applications [13][14][15][16][17][18].Remote sensing instruments are currently a widely used technology for analyzing LULC change in diferent regions of the world [19][20][21].Landsat, Sentinel-2, MODIS, the main product of wide coverage, and high-repeatability remote sensing imagery were recently used for LULC and water resource applications [7,22,23].
NDVI (Normalized Diference Vegetation Index), the most common Vegetation Index among remote sensing techniques, can indicate the growth status, type, and biomass of vegetation and has a linear correlation with vegetation cover.NDVI is a vegetation indicator that is considered the best indicator of vegetation growth and cover [24][25][26].Currently, scientists have obtained NDVI data from various satellites such as Landsat, Sentinel-2, and MODIS.Vegetation Cover Index analyzes were not conducted in the study area showing changes in NDVI in the area, despite some LULC studies showing slight LULC changes.
Te remote-sensing images used in this article are from the USGS Earth Explorer and NASA websites.Te Moderate Resolution Imaging Spectroradiometer (MODIS) is a satellite that collects remote sensing data used by scientists to monitor, model, and assess the efects of natural processes and human actions on the Earth's surface.MODIS has two sensors orbiting the Earth: the Terra (EOS AM) satellite, launched by NASA in 1999, and the Aqua (EOS PM) satellite, launched in 2002 [27][28][29].
Change in LULC is the major challenge in Ethiopia's hydroecological studies particularly in the study region [30].Te anthropogenic factor is the main source and results in a signifcant change in the LULC in Ethiopian watersheds.Key anthropogenic activities, including tree felling, land conversion to agriculture, and human settlement, are the major causes of LULC change.Tese changes in LULC were also recently studied by several researchers in Ethiopia [31][32][33][34].Some researchers tried to study the LULC change in the study region [30,35,36].However, the spatial scope of these works is limited and the period is also very short.Also, the methods used in previous studies do not address fnished products for remote sensing applications such as MODIS.None of the researchers have yet dealt with the LULC change of western catchments considering the spatiotemporal patterns of LULC.Te region considered in this particular study is an area ecologically rich and sensitive to climate and LULC changes.Te novelty of this study consists in comparing the methods (LULC and NDVI techniques) and their spatial and temporal coverage of the study region.In addition, the study updates LULC and vegetation information by mapping changes in magnitude, trend, and spatial distribution.Identifying and mapping the LULC change of a given area are critical for understanding the vulnerable area creating a sustainable ecosystem for hydroclimatic management.Terefore, this study aimed to analyze the trends and changes in LULC and vegetation dynamics in western Ethiopia with MODIS satellite imagery between 2001 and 2020.Te results of this study could provide frsthand information on the trend, extent and magnitude of LULC, and Vegetation Index in the study region.Furthermore, the study provides foundational information for LULC and vegetation-related information for regional and national governments, policymakers, stakeholders, and the local community.

Study Area Description.
Te study area is located in western Ethiopia and includes the Gambela Regional State in the west and some areas in the Oromia Regional State in the east.Te area is part of the Baro-Akobo River Basin, which is one of the tributaries of the Nile River system.Te catchment is located at latitudes 6 °50′ and 8 °43′ north and longitudes 33 °00′ and 35 °51′ east and covers an area of over 40889 km 2 .Te area receives the highest amount of mean annual precipitation (1600 mm) and humidity is the climatic character of the region.Te relative humidity in the region ranges from 62.8 to 77.9% while the mean temperature ranges from 18.08 to 28.8 °C.Te area is characterized by more than fve rainy months between March and October.Forest and savanna are the most important LULCs dominating more than 70% of the area.

Data Acquisitions.
Two fles were obtained and evaluated for this study: MCD12Q1 Land Cover Type (Collection 5) and MOD13Q1 Vegetation Indices.MODIS products were accessed from freely available data from the Earth Observing System Data and Information System (EOSDIS) Internet page.Te MODIS Land Cover Type product includes many categorization systems that describe land cover attributes generated from observations spanning a year's worth of Terra and Aqua data input.We used the most current Collection 5 MODIS Global Land Cover Type product, which was released in 2020 and spans the years 2001-2020 with a spatial resolution of 500 m.For the years 2001-2022, the MODIS 16-day composite vegetation index product was examined.It has a spatial resolution of 250 m and comprises the Enhanced Vegetation Index (EVI) and the Normalized Diference Vegetation Index (NDVI).Te overall method of procedure of the study is presented in Figure 1.
LULC data (Figure 2), NDVI data (Figure 3), climate data (Figures 4(a) and 4(b)), and topographic data (Figure 4(c)) were used for this particular study.Te Moderate Resolution Imaging Spectroradiometer (MODIS) product was chosen to analyze the spatial and temporal trends of LULC and NDVI.Te MODIS LULC data from 2000 to 2020 were downloaded from the United States Geological Survey (USGS).For NDVI data, the MOD13Q1 MODIS was used from the NASA website (Table 1).MODIS satellite was launched and developed by NASA in 1999.MODIS is a key instrument onboard the Terra (originally known as EOS AM-1) and Aqua (originally known as EOS PM-1) satellites [27,28].For 2001, present, the MODIS Land Cover Type product (MCD12Q1) provides worldwide maps of land cover at yearly time steps and 500-m spatial resolution.Te new data from the improved Collection 6 MODIS Global Land Cover Type product were utilized in this work to increase and correct spatial resolution utilizing the EVI, LST, and NBAR measures [37].

Data Processing.
Te data products obtained from MODIS observations describe aspects of the land, seas, and atmosphere that may be utilized for local to global scale investigations of processes and trends.MODIS NDVI data from the MOD13Q1 product ofer a Vegetation Index (VI) value per pixel.Te Normalized Diference Vegetation Index (NDVI), also known as the Continuity Index to the existing NDVI, is generated by the National Oceanic and Atmospheric Administration-Advanced Very High-Resolution Radiometer (NOAA-AVHRR).Te Scientifc World Journal We frst collected the available data of LULC and NDVI products of MODIS from the website.Te analysis in this study was carried out over a period of two decades, commencing in 2001, and used LULC and NDVI data to measure changes in LULC and vegetation change typically over a period of roughly 5 years.MODIS LULC data were utilized from 2001 to 2020 based on availability, whereas MODIS NDVI data were used from 2001 to 2022.All the downloaded MODIS pictures were collected, preprocessed, mosaicked, projected, and processed using Arc GIS tools.MODIS13Q1 datasets with the same temporal resolution (16 days) and spatial resolution (250 m) were obtained from the MODIS website from 2000 to 2022.Te time series data were checked for noise and parameter reconstruction to obtain correct parameter consistency in each pixel of the datasets for analysis.Te previous land cover data of the study region were collected from the Ethiopian Geospatial Agency.
Te processing stage was carried out using Arc GIS tools after getting the data from the website.Te frst stage was to ensure that the study's data covered the appropriate area.Mosaicking of separate pixels for each year of the dataset is the next step of the process.Te layers were trimmed to the research region after mosaicking to a new raster dataset.Reclassifying the image was the next step of the analysis which helped to reduce the number of classes.Tere were 17 layers in the MODIS LULC image dataset.Among the 17 layers, we have only 10 layers available in the study region.Table 2 shows the outcomes of marginal classes in the study region.Te steps of mosaicking, projecting coordinating, and classifying were done for each year from 2001 to 2020.Mapping the LULC was the next step of the analysis (Figure 1).
Te changes and dynamics of diferent land cover classes were investigated using net change analysis and crosstabulating matrix analysis methods.Investigating the net

Accuracy Assessment of LULC.
Te comparison between the classifed image and the ground truth data can be verifed by accuracy assessment.Accuracy assessment of a classifed image is important for analyzing LULC changes and determining the acceptability of the classifcation process [39][40][41].To check the accuracy of the LULC interpretation, 82 points were randomly selected in the study area and verifed using Google Earth software.Te initial step was to generate randomly generated points from the categorized LULC picture.For analytical assessment, these points are used as a reference between the categorized picture and the ground truth.Te producer's accuracy (PA) (equation ( 3)), user's accuracy (UA) (equation ( 2)), overall accuracy (OA) (equation ( 1)), and kappa coefcient (K) (equation ( 4)) are the parameters used to evaluate the LULC classifcation.Te equations for each assessment are presented as overall accuracy � total no. of correct classified pixels No reference pixel * 100, user acuracy � no. of correct classified pixels in each class Total no. of classified pixels in that class * 100, procedure acuracy � no. of correct classified pixels in each class Total no. of referenced pixels in that class * 100, kappa coefficient (T) � (total sample * total corrected sample) − (column total * row total) total sample 2 − (column total * row total) * 100. (4)

Change Detection.
Te change detection is the process of fnding diferences in the state of an object by observing it at diferent points in time [21,42].Te goal of change detection is to examine the variability in LULC recorded during a certain time period connected to a specifc location.Te general steps of the change detection model were as follows: data collection, preprocessing of remotely sensed data, image classifcation using appropriate logic, change detection using postclassifcation analysis, and vectorization of the raster change layers so that they could be manipulated in GIS tools.Te change from one LULC class to another class between 2001 and 2020 was quantifed.Te distribution of the changes that occurred for the specifed period is also presented in the matrix tables.A change matrix plots LULC changes from 2001 to 2020, assessing overall changes in LULC classes [4,9].Te examination of alterations in LULC from 2001 to 2020 was also demonstrated through the measurement of variation, the amount of change, and the distribution of maps that highlight the spatial and temporal diferences.To investigate the maximal range of change, spatial and temporal LULC alterations were detected for LULC types.Temporal changes were investigated using the total diference in areas for each LULC between two specifc times, whilst spatial changes were investigated on a pixel-bypixel scale to assess intraclass variations throughout these periods.Te pattern of shifts was also exhibited to clarify the analysis of changes in LULC and NDVI.Te change in the area of each LULC identifed in Arc GIS was calculated using Excel.Te categorized picture is manipulated using Arc GIS software for the frst phase of

Results
Arable land, evergreen broadleaf forest, deciduous broadleaf forest, grassland, permanent wetland, savanna, woody savanna, mixed forest, urban lands, and water bodies are the main LULC classes for the catchment (Table 2).Savanna and evergreen broadleaf forests are the main LULCs covering about 75% of the area.However, built-up areas and water bodies are the LULCs with the lowest area coverage (Figure 2).

Spatial Distribution of LULC and NDVI.
Te catchment is dominated by savanna grasslands as the main LULC type (Figure 2).Savana dominated most of the western plains of the catchment for the past two decades.Evergreen broadleaf forests are the other LULC type that represents the southeastern portions.Distribution of the deciduous broadleaf forest is along the western portion of the evergreen broadleaf forest in the southeast.Cropland is distributed in an uneven and pitted pattern in the central part.Urban and built-up areas, mixed forests, and water bodies are the LULCs with the lowest area coverage.Forest areas are areas with higher elevations and the lowlands are dominated by savanna.Te verdant regions exhibit a superior NDVI score of over 0.3, while the zones ranging from yellow to red depict lower NDVI values of under 0.3 (Figure 4).Typically, the western zone of the watershed encompasses a reduced NDVI value when compared to the southeastern forest areas, which boasts an NDVI value almost reaching 1.0.(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008).Te most important change between 2001 and 2008 was the change in large-scale savanna LULC to grassland (1026 km 2 ), followed by the change from evergreen forest to woody savanna (623.2 km 2 ) and grassland to savanna (311.4 km 2 ) (Figure 5).In general, savanna, grassland, and cropland are the main LULCs where changes occurred in a shift to each other.Te reduction and increment of each LULC in the given period are presented in Table 3. (2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015).Te changes from largescale savanna to grassland (962.5 km 2 ) and from evergreen forest to woody savanna (545.7 km 2 ) are the largest changes as in the previous trend (Figure 6).Savanna to deciduous forest (395.1 km 2 ), grassland to savanna (319.9 km 2 ), and savanna to woody savanna (223.1 km 2 ) are the other major changes covering large areas.Tere are also signifcant changes such as woody savanna to grassland (146.3 km 2 ), grassland to woody savanna (142.5 km 2 ), woody savanna to evergreen forest (134.7 km 2 ), and deciduous forest to savanna (132 km 2 ).Table 3 shows a decrease of 1035.5 sq•km in the savanna area and a decrease of 468.7 sq•km in the evergreen forests.Grasslands and forested savannas were the areas with the greatest increase between 2008 and 2015.

LULC Change (2015-2020).
Here, the switches from grassland to savanna (878.8 km 2 ) and from savanna to grassland (587.8 km) are the most important changes.Te change from the evergreen forest (350.1 km 2 ), grassland (341.2 km 2 ), and savanna (313.4 km 2 ) to woody savanna is the main feature of the LULC change in this particular period (Figure 7).In addition, there was a change from woody savanna to evergreen forest (314.5 sq km).Te change in coverage area between 2015 and 2020 is presented in Table 4. Tere have been major changes to the savannas and woody savannas LULC classes.Change (2001Change ( -2020)).Te detection of changes revealed that the savanna has changed into extensive grassland (1739 km 2 ) and deciduous forest (1640 km 2 ) (Figure 8).Te other major change is the change of EBF to a woody savanna (1110 km 2 ).Savanna to woody savanna, grassland to woody savanna, deciduous forest to evergreen forest, cropland to savanna, and savanna to evergreen forest are the other major changes at 513, 454, 407, 369, and 206 square kilometers, respectively.Te rate and percentage change between 2001 and 2020 are shown in Table 4).

LULC
Te changes in cover of the savanna, evergreen broadleaf forest, and croplands are the largest changes that occurred in the catchment, which decreased to −2975, −571, and −516 sq•km respectively, (Figures 9 and Table 4) between 2001 and 2020.In recent decades, there has been an increase in coverage for woody savannas, deciduous forests, and grasslands.Tere is no signifcant change for permanent wetlands, water bodies, and urban areas.
Savanna and forest are the major LULCs that cover a total of more than 70% of the area.Savannas are distributed in large areas in the west while forest dominated the southeastern part for the whole period of the analysis (2001 to 2020) (Figure 10).Te increment grasslands in the savanna-dominated areas of the northwestern coroner give 6 Te Scientifc World Journal a new appearance to the LULC map.Te sparse distribution of croplands in the savanna areas of the 2001 map also got a very tiny distribution in the 2020 map (Figure 2).

Trend of LULC (2001-2020).
Uniform trends of increase, decrease, and uniformity (no change) are cases where trend detection results were displayed.In other words, no inconsistent status is seen in the trend analysis of the catchment.Te LULC trend (Figure 11) shows that evergreen broadleaf forests, savannas, and cropland show a declining trend while deciduous broadleaf forests, woodland savannas, and grasslands are increasing from 2001 to 2020.Tere is no discernible trend in permanent wetlands, urban areas, and water bodies (Figure 11).Te changing trend analysis shows that perfect positive trends were seen for deciduous forests and woody savannas (Figures 11 and 12).Besides the change between 2015 and 2020, grassland shows a positive trend for most of the periods.For savannas and croplands, there was a negative trend in the past.With the exception of the most recent years (2015 and 2020), when they exhibit an upward trend, evergreen deciduous forests have historically shown a declining trend.

Accuracy
Assessment.An accuracy assessment was performed between the MODIS LULC data and Google Pro data from 82 random points taken from the study site.
Accuracy is checked using the most commonly used assessment method called kappa accuracy.Te matrix table (Table 5) displays the outcomes of precision verifcation by utilizing a randomly selected point from the LULC map and the Google Earth data.Te details of the method and user accuracy for every LULC parameter are presented in Table 6.

NDVI Characteristics.
Te NDVI values are indices that indicate the health of vegetation in a specifc area (Figure 3).Te Vegetation Index value generally increased from 2001 to 2022.However, the increase in value is not uniform, and it is diferent in dry and wet period trends (Figures 13 and 14).In the dry period (January), it increased between 2001 and 2008, decreased between 2008 and 2015, and increased again between 2015 and 2022 (Figure 13).Te wet period average value of NDVI shows an opposite trend to the dry period which decreased from 2001 to 2008, increased from 2008 to 2015, and again decreased from 2015 to 2022 with only a low range compared to dry period cases (Figure 14).Te classifed value of NDVI showed a similar general trend for dense vegetation area coverage.From 2001 to 2022, the vegetation cover and growth were improved by covering a large area.In addition to increasing the density of green vegetation in the southeastern part, the western corner areas were covered with green vegetation in late 2022 (Figure 3).Te Scientifc World Journal Te dry period NDVI shows that dense vegetation coverage increased from 27% 44% while barren land decreased from 35% to 10% between 2001 and 2020.During dry periods, a declining trend in the coverage of NDVI values (0.015 to 0.181) was observed between 2001 and 2022.On the other hand, during the same time period, the coverage area showed comparable growth for NDVI values greater than 0.27.
Te dry period analysis of NDVI shows that the value range of the NDVI_2001 value (0.0151 to 0.18) covers the largest area followed by the value of the NDVI value (0.36 to 1.0) (Figure 13).In the wet period, the NDVI value of 2022 shows the largest value covered by the NDVI value (0.361-1.0) followed by the NDVI value (0.27-0.36   ).Te wet period analysis of NDVI shows that the majority of the area is covered by a value of 0.36 to 1.0 for all the four periods (Figure 14).

Discussion
Land-cover maps are critical variables in studying envidynamics.Change in LULC is unavoidable since it is a result of social and economic progress.However, it frequently comes with a high environmental cost.Te expansion of agricultural land, for example, causes water pollution, sedimentation, erosion, and loss of biodiversity.Changes in land cover are usually caused by human activity such as urbanization, agricultural expansion, and deforestation, as well as natural occurrences such as wildfres, foods, and desertifcation [43].Changes in land use and land cover (LULC) have recently been one of the most important and persistent causes generating changes in the Earth's land.Tis study articulates the LULC and NDVI changes for the last two decades using MODIS satellite images.Te study discovered that LULC classes, namely, evergreen broadleaf forest, savanna, and water bodies, have reduced at the expense of natural vegetation covers such as shrubland and woodland.Change detection analysis was performed in order to understand the changes for the last period of time.Change detection analysis was performed in the Arc GIS environment to understand how a given area has changed between diferent time periods.We studied MODIS land cover categorization products from 2001 to 2020 and compared them to the produced MODIS NDVI trends to determine which land cover classes show trends and hence may be infuenced by the change.
Net change analysis and cross-tabulating matrix analysis methods were used to analyze the changes and dynamics of various land cover classes.Investigating the spatial distribution maps of diferent land cover classes and establishing their proportion of the total area of the research, it was discovered that the savanna (47.8%) and evergreen broadleaf forests (27.7%) are the dominant LULCs encompassing around 75% of the area followed by grassland, woody savannas, and deciduous broadleaf forest with a coverage of 8.7% 8.2%, and 7.3%, respectively.However, cropland, permanent wetlands, mixed forests, built-up regions, and water have the least amount of coverage of the total area.Te trend of LULC in the research region difers by LULC type from 2001 to 2020.Deciduous broadleaf forests, mixed forests, woody savannas, grasslands, and permanent wetlands all showed an increasing tendency.However, there was a declining trend in evergreen broadleaf forests, savannas, water bodies, and croplands.
Te dynamics of transition from one land cover type to another were studied, revealing a reciprocal link between various LULC types.Te rise in one LULC area was connected with a reduction in another LULC type.Decreased evergreen broadleaf forest, for example, is connected with increased deciduous broadleaf and mixed forest land.Te decrease in savanna coverage may be related to an increase in grassland area.Tese reciprocal relationships were also observed in water bodies and permanent wetlands; a decrease in water body coverage may be associated with an increase in wetland regions.
In all four independent eras of investigation, the signifcant change in LULC was caused by a shift from the savanna to grassland.From 2001 to 2020, the transition from  Te Scientifc World Journal savanna to grasslands in each epoch diminishes the area covered.However, the reciprocal change in grassland to savanna LULC increased with a narrower range from 2001 to 2020 than the savanna to grassland shift.Tis indicates that large areas of have changed to grassland nature of LULC over the past decades.Te other major change was the change from evergreen broadleaf forest to woody savanna which increases from time to time.Te fndings also demonstrate that the evergreen forest of south-western Ethiopia changes its cover from time to time, posing a severe environmental risk to the region and the country as a whole.In the research area, the greenness of vegetation increases with precipitation.Te range of NDVI values for the dry period is reasonably spread.During the rainy season, however, most of the land turns green, with an NDVI value greater than 0.36.Te greenness also follows the gradient of the slope.High-elevation places ofer more lush green fora than gradual sloping lowlands.
Te Vegetation Index results show that larger NDVI values occurred in forest and green vegetation areas, and the lowest value was found in barren and built-up uplands.Te evergreen, deciduous forest, and mixed forest areas of the southeastern parts had NDVI values ranging from 0.30 to 1.0.Most of the savanna and woody savanna areas have NDVI values of 0.10-0.40.Wetlands and the water bodies have NDVI values of −0.05 to 0.05.Barren and built-up lands have NDVI values in the range of 0.050-0.14.In addition to increasing the density of green vegetation in the southeastern part, the western corner areas were covered with green vegetation in late 2022 (Figure 3).Tis change in NDVI corresponds to the change in LULC at the end of 2020 which shows an increment of grassland coverage in the western part of the study.Te Scientifc World Journal 13 Tere are a number of similar studies that show comparative results of changes in LULC in Ethiopia as well as abroad.A study by Moisa et al. [33] on the Anger River Subbasin, western Ethiopia, shows results corresponding to our study.Other similar studies also have a corresponding result in western Ethiopia [43][44][45].Tere are several studies on the database that have principles analogous to the present research.Te work is similarly analogous to Mahmoudi et al.'s [38] study on detecting land cover changes in Baluchistan using the MODIS Land Cover product.Another comparable study was conducted by Somayajula et al. [46] on classifcation and validation of spatio-temporal changes  14 Te Scientifc World Journal in land use/land cover and land surface temperature of multitemporal images.
To manage overall environmental resources, an understanding of LULC variability in a given location is required.Te suggested technique has numerous applications in land resource management.Tis sort of research provides frsthand knowledge regarding the surface covering.Te LULC variability maps provide information for policymakers and general river basin planning in an area.Te diference between this research and others is related to its areal coverage as well as its spatial and temporal resolution.

Conclusions
Te present study shows how LULC and vegetation variability have changed signifcantly in western Ethiopia over the past two decades.Tis study evaluated and checked changes in LULC and NDVI patterns in western Ethiopia using MODIS satellite data from 2001 to 2020.Using MCD12Q1 and MOD13Q1 from 2001 to 2020 as data sources, this study examined the annual diferences in LULC and NDVI for diferent years.Te overall classifcation accuracy of LULC for the classifed image from 2001 to 2020 was found to be 85.4% and the overall kappa statistics was 81.2%.Te result of this study reveals that the major LULC in the catchment is savanna and forest.Te cropland, savanna, and water bodies were reduced by 90%, 13%, and 19%, and woody savanna, deciduous forest, grassland, and permanent wetland increased by 119%, 52%, 45%, and 37%, respectively, over the past decades.Te results also show that some water body areas have changed to permanent wetlands.Te trend result showed that deciduous forests, woody savannas, and grasslands has increasing trends through 2001 to 2020.But the evergreen deciduous forest shows a decreasing trend for the given period.Te most notable change was the change from savanna to other LULC classes such as grassland and forest.Te result also showed that the NDVI value of the catchment was −0.046 to 0.98 in the driest month and it ranges from −0.093 to 0.99 for the wettest month.In the wettest month, 80% of the area has an NDVI value range of 0.36 to 1.0.Tis shows that 80 of the area is covered by green vegetation for the rainy season.Te researcher concludes his research by stating that the catchment has achieved dual coverage of LULC and is losing that coverage that needs to be preserved for the future.Tus, this investigation foresees that the results could furnish insights into those responsible for overseeing water and land, as well as policymakers, in order to ensure the enduring management and growth of the area's natural assets.
Detection.Here, we describe the changes between classes (conversion from one LULC class to another) and within classes.In this study, LULC change is detected as a change in LULC labels between 2001 and 2020 MODIS imagery.LULC change detection for the area is presented in the following tables and fgures.Te results are presented in a change matrix table that contains important information about the change from one LULC class to another class.Te change matrix showing the land cover changes between 2001 and 2020 was generated from classifed images of the respective periods and used to assess the overall changes in the LULC classes between 2001 and 2020.

Figure 4 :
Figure 4: Te climatic and topographic map of the study area.

Table 1 :
Characteristics of MODIS LULC and NDVI data.
[38] it does not provide precise knowledge about the conversion of one class to another.In other words, determining which land cover class was changed to another or which class caused fundamental changes in other land cover classes is impossible[38].

Table 2 :
Description of LULC of the catchment.Change detection analysis involves converting raster images into shapefles, merging several layers, and calculating the area of each LULC.Tese steps were performed for each analysis period(2001, 2008,  2015, and 2020).Te next step was to perform change detection using the intersection tool of Arc GIS.
).Generally, the greenness of vegetation increases from 2001 to

Table 6 :
Procedure and user accuracy assessment results of LULC 2020.